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ABSTRACT 

By analogy with the minimum-mass solar nebula, we construct a surface- 
density profile using the orbits of the 26 precise-Doppler planets found in 
multiple-planet systems: S = 2200 gcm~ 2 (a/1 AU)~^, where a is the circum- 
stellar radius and j3 = 2.0 ± 0.5. The minimum-mass solar nebula {(3 = 1.5) is 
consistent with this model, but the uniform-a accretion disk model (/3 ~ 1) is 
not. In a nebula with (3 < 2, the center of the disk is the likely cradle of planet 
formation. 

Subject headings: planetary systems: protoplanetary disks — solar system for- 
mation 

1. INTRODUCTION 



Our understanding of the mass profile of the disk from which the planets of the solar 
system formed begins with the minimum mass solar nebula (MMSN). Imagine adding just 
enough light elements to each planet so that all the planets reach solar composition, then 
smearing the augmented planets into contiguous nested rings (Edgeworth 1949; Kuiper 1956; 
Safronnov 1967; Alfven & Arrhenius 1970; Kusaka et al. 1970; Lecar & Franklin 1973). Fit- 
ting a power-law to the resulting surface-density distribution leads to the well-known MMSN 
prescription, £ = E (a/1AU) -3 / 2 , where £ is the surface density and a is the circumstellar 
radius (Weidenschilling 1977; Hayashi 1981). 
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The estimates of Weidenschilling (1977) for the masses of the planets augmented to 
solar composition suggest that S ~ 4200 gem -2 . Using the same table of solar abundances 
(Cameron 1973), but making some additional assumptions about the importance of the snow 
line, Hayashi (1981) derived S = 1700 gem -2 . Since the snow line's location and role in 
giant planet formation are uncertain (Sasselov & Lecar 2000), we prefer the above formula 
based on Weidenschilling (1977) for the purposes of this paper. 

Elaborate calculations of planet formation (e.g. Trilling et al. 1998; Ida & Lin 2004) 
regularly invoke the MMSN or a E ~ a -1 surface density distribution based on uniform-a 
accretion disk models (Shakura & Sunyaev 1973). Such calculations can depend strongly 
on the assumed surface-density distribution. For example, the isolation mass produced 
by runaway planetesimal growth goes as £ 3 / 2 (Lissauer 1987). The locations of secular 
resonances in the nebula are especially sensitive to the surface-density distribution (Lecar & 
Franklin 1997; Nagasawa & Ida 2000). 

Now, newly discovered extrasolar planets outnumber the solar system planets. This pa- 
per attempts to update our picture of the minimum-mass nebula taking the extrasolar planets 
into account; we derive a surface- density profile analogous to the MMSN using the orbits of 
the 26 extrasolar planets found in multiple-planet systems. Perhaps this new construction — a 
minimum-mass extrasolar nebula (MMEN) — can serve future studies of planet formation. 

2. DATA 

The extrasolar planets discovered so far by precise-Doppler methods contain nine two- 
planet systems and two three-planet systems. We obtained the orbital data for these planets 
from the website http://exoplanets.org; see the review by Marcy & Butler (2000) for back- 
ground on radial-velocity planet-searches. Table 1 summarizes the relevant data on these 
systems including measured stellar metallicities from Santos et al. (2004). Of the eleven 
listed extrasolar planet host stars, eight are G type. 

Besides multiple systems inferred purely from radial velocity observations, the table also 
includes e Eridani. This system contains one planet inferred from radial-velocity measure- 
ments to have a semimajor axis of 3.3 AU and a second planet inferred from the structure 
of a circumstellar debris disk imaged in the submillimeter (Greaves et al. 1998; Quillen & 
Thorndike 2002; Kuchner & Holman 2003). Quillen & Thorndike (2002), estimate that this 
second planet has a semimajor axis of 40 AU and a mass of 10~ 4 times the mass of e Eridani 
(0.8 M ). Since the existence of these planets is still debated, we perform our analysis both 
with and without them in the sample. 
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Table 1: Multiple Planet Systems 



Star 


Spectral 


[Fe/H] 


Semimajor Axes 


So 


P 




Type 




(AU) 


(gem J 




Sun 


G2 V 


0.0 


5.20, 9.54, 19.19, 30.07 


4225 


1.78 


55 Cancri 


G8 V 


0.33 ±0.07 


0.115, 0.24, 5.9 


739 


2.42 


v Andromeda 


F8 V 


0.13 ± 0.08 


0.059, 0.829, 2.53 


2670 


1.50 


GJ876 


M4 


0.0 a 


0.13, 0.21 


8634 


2.00 


HD 38529 


G4 IV 


0.40 ±0.06 


0.129, 3.71 


705 


1.58 


HD 82943 


GO 


0.30 


0.728, 1.16 


8907 


2.00 


HD 169830 


F8 V 


0.21 ±0.05 


0.81, 3.60 


2815 


2.00 


HD 12661 


G6 V 


0.36 ±0.05 


0.82, 2.6 


3013 


2.00 


HD 168443 


G5 IV 


0.06 ±0.05 


0.295, 2.87 


3374 


1.65 


HD 74156 


GO 


0.16 ±0.05 


0.294, 3.40 


1012 


1.63 


47 Ursa Majoris 


GO V 


0.06 ±0.03 


2.89, 3.73 


7026 


2.00 


HD 37124 


G4 V 


-0.38 ±0.04 


0.54, 2.5 


1949 


2.00 


e Eridani 


K7 V 


-0.13 ±0.04 


3.3, 40.0 


715 


2.00 



a From http:/exoplanets.org 

Radial velocity surveys are generally biased against detecting low mass planets at large 
periods. Precise-Doppler surveys can presently detect radial velocity variations with ampli- 
tude > 3ms± In terms of the planet's orbital period, P, and the ratio of the mass of the 
planet to the mass of the star, /i, this detection limit corresponds to 

/ /sin*>io- 4 (— y /3 fj^y 1/3 (i) 



v lyr; \M &/ 

where % is the inclination of the normal to the planet's orbital plane to the line of sight. 

Figure 1 shows the orbital period and fisini of each planet. Open circles indicate the 
planets in two-planet systems. Filled circles indicate planets in three-planet systems. The 
triangles indicate the e Eridani planets. The solid line shows the 3 ms -1 detection limit, 
calculated for a 1 M star. For the purpose of Figure 1, we assume sini for e Eridani is 0.5 
based on the inclination of the disk. 

This figure suggests that the detection bias in this data set is likely to be small; all of 
the precise-Doppler planets in our data set could have been detected out to the period limit 
(roughly 10 years), except for 55 Cancri c, which could have been detected only out to a 
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Fig. 1. — Orbital period and /xsini for the extrasolar planets in multiple-planet systems; 
is the ratio of the mass of the planet to the mass of the star. The solid line shows the 3 
ms _1 radial- velocity detection limit for a 1 M star, comfortably below yusini for most of 
our sample. 
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period of 6 years. However, other biases might lurk in the process of reporting and publishing 
the exoplanet data, and naturally, detecting planets via disk perturbations suffers from a 
special set of biases. Our analysis bears repeating when more data become available. 

3. ANALYSIS 

Constructing a MMEN from these data requires associating a surface density with each 
extrasolar planet. To derive these surface densities, we first must estimate an augmented 
mass for each planet — the mass each planet would have if it accreted just enough additional 
gas with the right chemistry so that it attained solar composition. We can convert the 
augmented mass to a surface density based on the spacings of the planets in multiple-planet 
systems. 

Jupiter has a mass of 318 and roughly 2-4 times solar abundances of C, N and 
S (Gautier et al. 2001), suggesting that its augmented mass should be roughly 1000 Mm. 
Estimates for the augmented masses of Saturn, Uranus and Neptune are all consistent with 
this value, even though their actual masses span a factor of ~ 20 (Weidenschilling 1977). 
This trend supports a planet formation scenario in which the observed range of planet masses 
arises mainly from variations in volatile depletion. 

We chose 1000 for the augmented masses of most of the extrasolar planets. If 
we try to construct the minimum mass solar nebular from the orbits of Jupiter, Saturn, 
Uranus and Neptune under this approximation, we get £ = 4230 gcm~ 2 (a/1 AU)~ L8 , in 
rough agreement with the MMSN of Weidenschilling (1977). A few planet candidates have 
measured minimum masses greater than 1000 M^: HD 168443b (msini = 7.73 Mj), HD 
168443c (msin? = 17.2 Mj), HD 74156c (msini = 8.21 Mj), v Andromeda d (msini = 
3.75 Mj), and 55 Cancri d (msini = 4.05 Mj). We used the measured minimum masses, 
msini, as the augmented masses for these planets. Using msini for all the augmented 
masses would substantially increase the dispersion in the surface densities. 

To derive surface densities, we divided each planet's augmented mass by the area of an 
appropriate annulus. Here we harnessed the multiple planet systems; we used the separation 
of the planets in log semimajor axis space as the full width of each annulus. In other words, 
we chose the boundaries of the annuli to be the geometric means of the adjacent semimajor 
axes. Weidenschilling (1977) found that using an arithmetic mean or a geometric mean 
generated indistinguishable results for the solar system. However, the larger separations of 
the extrasolar planets seem to call for the use of the geometric mean, given the usual picture 
of feeding zones for planetary accretion that span widths measured in units of the local Hill 
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radius (e.g. Kokubo & Ida 1998). The extrasolar planets exhibit a much larger range of 
orbital eccentricities than the solar planets (Marcy & Butler 2000); we do not attempt to 
take this factor into account. 

Figure 2 shows the surface densities we associate with each planet. Open circles indicate 
the surface densities associated with the planets in two-planet systems. Filled circles indicate 
the surface densities associated with the planets in three-planet systems. Also plotted are 
surfaces densities associated with the e Eridani planets, marked by triangles. 

Table 1 reveals some of the subtleties of this procedure. It lists the parameters S and 
P for minimum-mass nebulae fitted separately to each system. The table shows that (3 = 2 
for any two-planet system where both planets have the same mass; this limitation of the 
algorithm prevents us from taking the separately-fitted nebulae too seriously or subdividing 
the data set in too much detail. 

The table also shows that the metallicities of the host stars range over 0.67 dex. At 
first, one might wonder how this spread in metallicity should affect the augmented masses. 
However, it turns out that for power-law models fit without any attempt to take account of 
stellar metallicity, So shows no correlation with stellar metallicity, only scatter. This lack 
of a discernible trend suggests that metallicity corrections are beneath the fidelity of our 
method. 

Table 2 shows S and (3 for power laws fit to the data taken as a whole and to 
subsets of the data. The power law best fit to the whole data set has the form E = 
2200 gcm~ 2 (a/1 AU)~ 2 - . This form for the MMEN has a total mass of 8 Jupiter masses in 
the region 0.3-30 AU, compared to 31 for the MMSN. This difference may result from the 
present inability of radial velocity techniques to measure orbital periods longer than ~ 10 
years; we may have only found 1/3 of the mass we will eventually find. 

Figure 2 compares the above form for the MMEN (solid line) with the standard MMSN, 
£ = 4200 gem -2 (a/1 AU)~ 3 ^ 2 (dashed line). Assuming the logarithms of the surface densi- 
ties have uncertainties of ±0.3, the uncertainty in the exponent of the power-law fit to the 
entire data set is (3 = 2 ± 0.5. Given these assumptions, the standard MMSN is consistent 
with the exoplanet data. 
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Fig. 2. — Extrasolar planet surface densities. Filled circles are three-planet systems and 
solid circles are two-planet systems. Two triangles indicate the e Eridani system. The best 
fit power law is a~ 2 (solid line), though a~ 3 / 2 (dashed line) is also marginally consistent with 
the data. 
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Table 2: Best fit surface-density profiles. 



Data Set 


S (g cm 2 ) 


P 


Two-Planet Systems Only 


2990 


1.87 


Three-Planet Systems Only 


739 


2.42 


G Star Exoplanets Only 


2150 


1.79 


All Exoplanets But e Eridani 


2500 


1.87 


All Exoplanets 


2220 


2.01 


Solar System Giants Only 


4230 


1.78 


Exoplanets + Solar System Giants 


2490 


1.87 



4. DISCUSSION 

4.1. Other Nebular Models 

By combining the data from several exoplanet systems, we have implicitly assumed that 
protoplanetary disks have a common surface- density distribution. Such a common profile 
might stem from a steady-state model of the protoplanetary accretion disk, like the uniform-a 
accretion disk model (Shakura & Sunyaev 1973). 

However, the uniform-a accretion disk model (/3 ~ 1) is not consistent with the MMEN 
power law, and perhaps we should not expect it to be. This primary justification for this 
model is that a uniform a is the simplest assumption one can make about anomalous disk 
viscosity. However, the reason for angular moment transport in disks is highly debated; 
models of disk turbulence driven by magneto-rotational instability (e.g. Fleming & Stone 
2002) suggest that the angular momentum transport is highly non-uniform, and may not 
even be characterized by a viscosity. 

An alternative steady-state disk model is a massive circumstellar disk bordering on grav- 
itational instability. Spiral density waves may form in these disks that tend to redistribute 
angular momentum to maintain Toomre's Q at ~ 1 throughout the disk (Laughlin 1994). 
Or else, planets could conceivably form via direct gravitational collapse of such disks (e.g. 
Boss 2001), depending on the thermodynamic assumptions (Pickett et al. 2000). For a disk 
with uniform Q, the surface density profile is determined by the square root of the local 
disk temperature, so a variety of disk thermal models would yield similar surface-density 
distributions. These disks are at least an order of magnitude more massive than the min- 
imum mass nebulae we are considering. For example, a Q — 1 disk where the midplane 
temperature is T = 150 K (a/1 AU)~ 3 / 7 (Chiang & Goldreich 1997) has surface density 
E = 88000 gcm- 2 (a/l AlfT L93 . 
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Of course, the real solar nebula may seldom have been near a steady state; T Tauri stars 
are highly variable, and at any given stellar age, one can find stars both with and without 
accretion disks in the same cluster (e.g. Haisch et al. 2001). One-dimensional models of 
protoplanetary disk accretion more sophisticated than the uniform-a model (Gammie 1996; 
Matsuyama et al. 2003) generate time-variable surface density profiles. Ironically, these 
models often use the MMSN as a starting point. In general, a non-steady-state picture of 
the protoplanetary disk stretches our interpretation of the MMSN and MMEN; all of the 
minimum mass need not have been available as free gas in the disk at the same time! 

4.2. A Steeper Surface-Density Distribution? 

Though taken as a whole, the extrasolar planet data agree with the MMSN, Figure 2 and 
Table 2 show that some extrasolar planetary systems themselves do not. The (5 exponents 
for the three-exoplanet systems and for the solar system differ substantially; the minimum- 
mass nebula constructed from the three-exoplanet systems has more mass concentrated in 
the center, like the uniform-Q model mentioned above. These three-exoplanet systems are 
the ones for which our algorithm works the best. 

Conventional wisdom holds that giant planets form at 2-5 AU from the star, perhaps 
at the snow line, and that the close-in extrasolar planets migrated inward from where they 
formed (Lin et al. 1996). According to this picture, one might might expect the MMEN to be 
more centrally peaked than the MMSN. Perhaps the reason for the steep power law among 
the three-exoplanet systems is that these planets are more widely separated in semimajor 
axis than planets in the solar system; with our prescription, increasing planet separations 
steepens the power law. The planet-disk interactions that cause the planets to migrate to 
small a could create the wider separations; massive planets can clear out the region between 
them in a protoplanetary disk via slowly damped spiral density waves (Bryden et al. 2000), 
preventing more planets for growing in these regions. 

5. INSIDE-OUT PLANET FORMATION 

However, looking at things another way, the steeper nebular surface density distributions 
consistent with the MMEN could have dramatic implications for where planets are formed. 
Let us briefly indulge this alternative viewpoint. 

An annulus of the protoplanetary disk with width Ar has mass m = 27r£rAr. If this 
annulus is a feeding zone with Ar oc r and if X oc r _/3 , then m oc r 2_/3 . For the minimum- 
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mass solar nebula {(3 = 1.5) and uniform-a disk model ((3 ~ 1), m increases with radius, 
suggesting that more mass is available for planet-formation at larger distances from the star. 
This tale — an exercise in circular reasoning — has been a standard argument for why giant 
planets ought not to form predominantly at the center of the disk. 

However, for (3 > 2, m decreases with radius, suggesting that the center of the disk is 
the most natural site for planet formation. Perhaps giant planets predominantly form in the 
center of their disks and migrate outwards. 

Indeed, suitable outward migration mechanisms for giant planets exist: Type III migra- 
tion (Masset & Papaloizou 2003) and planet-planet scattering (Thommes et al. 1999; Adams 
& Laughlin 2003). Type III migration begins when the planet grows larger than a threshold 
mass sufficient to begin opening a gap in the disk (Arymowicz 2004). Some planet-planet 
interactions may also be required to produce the observed range of extrasolar-planet orbital 
eccentricities. 

Outward Type III migration will halt at the outer boundary of the disk. Matsuyama 
et al. (2003) have suggested that the disk may develop an edge or a gap at ~ 5 AU when 
the photoevaporation time scale matches the viscous spreading time scale. The exact nature 
and location of this edge remain uncertain. However, the disk probably photoevaporates 
first outside ~ 5 AU, providing a natural stopping place for outward planet migration. 

This inside-out planet-formation scenario has several attractive features. While no di- 
rect observational evidence has yet been found for a special structure at the snow line in 
protoplanetary disks, abundant new evidence from near infrared spectroscopy (Hartmann 
et al. 1993; Muzerolle et al. 2003) and interferometry (Millan-Gabet et al. 2001; Colavita et 
al. 2003) implies that T Tauri disks have dramatic opacity jumps at ~ 0.1 AU. Models of 
layered accretion disks (Gammie 1996; Armitage et al. 2001) suggest that even while most of 
the disk has Q » 1, the center of the disk could become gravitationally unstable, fostering 
planet formation by disk instability. Moreover, formation of the planets in the center of the 
disk could potentially solve two problems in the core accretion picture: 

Type I Migration of Giant Planet Cores 

Conventional simulations of the semimajor axis distribution of the extrasolar planets 
(Trilling et al. 2002; Armitage et al. 2002; Ida & Lin 2004) neglect Type I migration (Ward 
1997) on the grounds that no planets could form at the snow line if this mechanism operated. 
Although present calculations of the migration of the small bodies do not accurately include 
torques added in the coorbital region (Bate et al. 2003), these torques are likely to dominate 
the net Lindblad torque, not compensate for it, exacerbating the problem. However, several 
mechanisms have been suggested that could stop migration of cores in the center of the 
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disk (Lin et al. 1996; Kuchner & Lecar 2002; Terquem 2003; Laughlin et al. 2003; Nelson & 
Papaloizou 2003). 

The Missing Short-Period Pile-up 

Tabachnik & Tremaine (2002) and Tremaine & Zakamska (2003) pointed out that the 
distribution of planet orbital periods shorter than ~ 200 days is consistent with a power 
law, and that scenarios where planets migrate inwards and stop migrating at a special small 
orbital period should produce a much larger pile-up of planets at short periods than observed. 
Gu et al. (2003) invoked tidal-driven Roche-lobe overflow to remove this exaggerated pile-up 
of planets in 3-day orbits expected from Monte-Carlo migration models. However, most of 
the stopping mechanisms cited above are not specific to precisely three-day periods, and since 
tidal forces decay as cr 6 , tide-induced Roche-lobe overflow can probably not help prevent 
pile-ups at say, a 6-day period. If planets mostly migrated outwards when they grew beyond 
a threshold mass as Type III accretion models suggest, this missing pile-up problem would 
not arise. 

Two effects make the very center of the disk (0.05-0.1 AU) problematic for core acretion: 
heating of the planet by tidal interaction with the star (Gu et al. 2003) and the large 
critical core radius for gas accretion at high nebular temperatures (Bodenheimer et al. 2000). 
However, the region at ~ 0.2 AU could be fertile ground for planet formation by either core 
accretion or disk instability. The disk temperature there is potentially low enough that the 
critial core radius shrinks to where it is consistent with the inferred core mass of Jupiter and 
Saturn (Papaloizou & Terquem 1999; Ikoma et al. 2001) particularly if this region is shadowed 
by a puffed-up inner wall of the disk (Natta et al. 2001; Dullemond et al. 2001). Conceivably, 
the observed extrasolar planets all could have formed at 0.1-0.2 AU and migrated outwards 
to their current orbits. The solar system giants might also have formed via this process, 
though this possibility raises the question of the origin of the abundant low-temperature 
condensates in the envelopes of these planets (e.g. Gautier et al. 2001), a question still 
inadequately answered by any planet-formation model. 



6. CONCLUSIONS 

By analogy with the minimum mass solar nebula, we constructed the surface density 
distribution for a minimum-mass extrasolar nebula based on the orbits of the extrasolar plan- 
ets in two- and three-planet systems: £ « 2200 gem -2 (a/1 AU) -2 . The standard MMSN is 
consistent with the exoplanet data taken as a whole, though the uniform-a accretion disk 
(Shakura & Sunyaev 1973) is not. 
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In general, though the MMSN is consistent with the MMEN, our analysis illustrates 
that the logic that led to the MMSN supports a wide range of disk surface-density profiles, 
none of which can serve as a complete picture. The MMSN was based on Laplace's concept 
of the solar nebula as a smooth disk that broke up into rings that condensed into planets. 
In contrast, true protoplanetary disks are likely actively accreting and time variable. 

Our analysis led us to consider the ramifications of alternative surface-density distri- 
butions. We pointed out that in the centrally-concentrated nebular model suggested by 
the three-exoplanet systems, the center of the disk is the preferred site for planet formation. 
The extrasolar giant planets — and possibly the giant planets of the solar-system — could have 
formed at 0.1 — 0.2 AU and migrated outwards to their present orbits via Type III migration 
and planet-planet scattering. 

Our simple model is not intended to replace the MMSN, but to illuminate its shortcom- 
ings and to expose our solar-centrism. We must not remain chained to /3 — 1.5; minimum 
mass can sometimes mean maximum prejudice. 

Thanks to Ed Thommes, Jack Lissauer, Norm Murray, Andrew Youdin, Mike Lecar, 
Pawel Artymowicz, Scott Tremaine, David Stevenson, and an anonymous referee for com- 
ments. 
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